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We generalize a Landauer-type formula, using a real(8)spin-space Green function technique, to treat 
spin-dependent transport in quantum-coherent conductors attached to two ferromagnetic contacts. 
The formalism is employed to study the properties of components of an exact zero-temperature 
conductance matrix G, as well as their mesoscopic fluctuations, describing injection and detection 
of a spin-polarized current in a two-dimensional system where electrons exhibit an interplay between 
Rashba spin-orbit (SO) coupling and phase-coherent propagation through a disordered medium. 
Strong Rashba coupling leads to a dramatic reduction of localization effects on the conductances 
and their fluctuations, whose features depend on the spin-polarization of injected electrons. In the 
limit of weak Rashba interaction antilocalization vanishes (i.e., the sum of the matrix elements of 
G is almost independent of the SO coupling), but the partial spin-resolved conductances can still 
be non-zero. Besides spin-resolved conductance fluctuations and antilocalization, unusual quantum 
interference effects are revealed in this system leading to a negative difference between the partial 
conductances for a parallel and an antiparallel orientation of the contact magnetization, in a range 
of disorder strengths and for a particular spin-polarization of incoming electron with respect to the 
direction of Rashba electric field. 

PACS numbers: 72.25.-b, 72.25.Rb., 73.23.-b, 85.75.Hh. 



I. INTRODUCTION 

Due to its fundamentally quantum nature, electron 
spin has very limited choice of interactions through which 
it can couple to the environment. Thus, its polarization 
can survive various forces in metals or semiconductors, 
while charge carriers undergo many scattering events, for 
long enough time to allow for envisioned quantum tech- 
nologies that manipulate spin, such as spintronics^'^ or 
solid-state quantum computing with spin-qubits.'^ Spin- 
tronic engineering of spin-polarized currents, combined 
with conventional electronics that manipulate electron 
charge, offers an exciting prospect to assemble infor- 
mation processing, storage, and communication on the 
same chip. Recent experimental advances have provided 
an impetus to study spin-polarized transport in two- 
dimensional electron systems, posing theoretical chal- 
lenges^ where the ubiquitous "factor of two" for spin de- 
generacy in usual transport formulas has to be replaced 
by a more involved analysis of the interplay between or- 
bital and spin degrees of freedom of the electrons. 

While quantum-coherence of spin is essential for spin- 
tronic qubits, as well as some of the proposed spintronic 
devices for classical information processing;^ modeling 
of the orbital kinetics of electrons that carry spin has 
mostly followed the traditional route of device engineer- 
ing in charge-electronics: the charge transport (ballistic 
or diffusive) is described by a semiclassical Boltzmann 
formalism treating electrons and holes as classical par- 
ticles where quantum effects enter only through their 
kinetic energy (being determined by the semiconductor 
band structure) Although many applications require 
devices operating at room temperature (where phase- 
coherence of orbital degrees of freedom is usually washed 
out), from the fundamental transport physics point of 
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FIG. 1: Graphical depiction of our lattice model for a generic 
two-probe quantum spintronic FM-Sm-FM device: two fer- 
romagnetic contacts (FM), that act as a spin injector (left) 
and a spin detector (right) are attached to a two-dimensional 
electron gas (e.g., in InAs), which is confined to the xy-plane 
by a Rashba electric field along the z-axis. The white ar- 
rows denote one possible magnetization configuration for the 
spin-resolved measurements (e.g., spin-up along j/-axis is in- 
jected and both spins are collected). Each site hosts a single 
s-orbital which hops to four (or fewer for surface atoms) near- 
est neighbors. The hopping matrix element is t (within the 
sample), ti, (within the FM leads), and tc (coupling of the 
sample to the leads). The leads are semi-infinite and con- 
nected smoothly at ±00 to macroscopic reservoirs biased by 
the chemical potential difference fiL — I^r = eV. 



view, it is intriguing to explore quantum interference ef- 
fects involving both the charge and spin of the electrons 
in nanoscale samples at very low temperatures < IK — 
a playground of mesoscopic physics over the past two 
decades.'' Moreover, better understanding of the role of 
conserved quantum coherence of both wave functions and 
spin states in the conductions process is crucial for the 
devices that involve spin-polarized resonant tunneling?^ 
Recent theoretical efforts have been exploring mesoscopic 
aspects of spintronics by extending standard techniques 



2 



(such as, the Landauer-Biittiker scattering approach to 
quantum transport of spinless particles in finite-size sys- 
temsSiili) to describe transport in the presence of spin- 
dependent interactions by taking into account quantum- 
coherent dynamics of both orbital and spin degrees of 
£reedom4*ii Even though it is possible to add local spin- 
dependent effects into the scattering matrix directly, such 
extensions become highly non-trivial in the most general 
case where spin effects occur all over the system, as in 
the presence of SO interactionsAi^ 

A particularly important interaction for semiconduc- 
tor spintronics is the Rashba type of spin-orbit cou- 
pling, ^'^ which is different from the more familiar impu- 
rity induced and position dependent one in metals. It 
arises from the asymmetry along the z-axis of the con- 
fining quantum well electric potential (generating a mag- 
netic field in the electron rest frame) that creates a two- 
dimensional (2D) electron gas (xy-plane) on a narrow- 
gap semiconductor surface. Since the Rashba coupling 
can be tuned by an external gate electrodcfi^ it is envis- 
aged as a tool to control the precession of the electron 
spin in the Datta-Das proposal for a field-effect spin tran- 
sistor.^ Thus, in such spin-interference devices intrinsic 
magnetic degrees of freedom can be controlled by elec- 
trical means. The spin-FET embodies a paradigm of a 
spintronic two-terminal device, such as the one shown 
in Fig. n (i) at the source (injector), information is 
stored as spins in a particular orientation (spin-up f 
or spin-down |); (ii) the spins, being attached to con- 
duction electrons, carry the information along a non- 
magnetic sample where they are manipulated through 
spin-dependent interactions; and (iii) the information is 
read off at the drain (detector). The successful realiza- 
tion of such devices demands resolution of some of the 
paramount fundamental problems in current spintronic 
research: all-electrical injection, at room temperature^ of 
fully spin-polarized currents into a semiconductori^ii^ (on 
the same footing as demonstrated in metallic devicesii 
where large currents of cold-electrons can be injected and 
detected via Ohmic contacts); efficient spin detection at 
a drain electrode (spin measurement is impeded by the 
same feature that generates long coherence times mak- 
ing it technologically interesting — weak coupling to the 
environment^^) ; and the engineering of spin-orbit cou- 
plings 

The paper commences with the discussion of differ- 
ent types of quantum coherence that can be encountered 
in spintronic structures (Sec. These intuitive no- 

tions are then rigorously elaborated within the tensor 
product formalism for the joint description of spin and 
orbital quantum states, thereby making it easy to in- 
clude effects of quantum coherence in the computation 
of transport quantities relevant for spintronics. The cen- 
tral such quantity is the conductance matrix G describ- 
ing spin-resolved transport measurements. We then in- 
troduces an efficient technique to compute G by devel- 
oping an implementation (based on lattice Green func- 
tionsMi formulated in the tensor product of real and 



spin space, rather than wave functions) of the Landauer- 
Biittiker formulation of linear response transport. Our 
formalism can accommodates both local and non-local 
effects of spin dynamics. We then reexamine in Sec. IIIII 
the salient quantum interference effects of mesoscopics — 
conductance fluctuations and (anti)localization — for the 
matrix elements of G, as well as their sums, cor- 
responding to different spin-resolved measurements on 
a ferromagnet-semiconductor-ferromagnet (FM-Sm-FM) 
system where disorder is gradually increased in the cen- 
tral non-magnetic region. This allows us to sweep 
through different transport regimes: from ballistic (where 
few scattering events take place within the sample of size 
L since mean free path satisfies £ > L) to strongly disor- 
dered (where mean free path loses its meaning). Recent 
quantum-coherent treatments of spin-polarized transport 
in two-terminal devices have been confined to clean sam- 
p^ggiijig because the Datta-Das original spin-FET pro- 
posal (which has ignited much of the interest for semi- 
conductor spintronics) requires strictly ballistic transport 
regime for the charge degrees of freedom. However, some 
level of scattering (from impurities, boundaries, or at the 
FM-Sm interface) is inevitable, and it is important to un- 
derstand properties of spin-polarized injection, conduc- 
tion, and degradation of well-defined spin-polarization in 
such environment.^*^ For example, peculiar quantum in- 
terference effects, arising from the entanglement of spin 
and orbital degrees of freedom^^ are revealed in our 
disordered FM-Sm-FM system when analyzing the rela- 
tionship between partial conductances (elements of G) as 
a function of the Rashba SO coupling, disorder strength, 
and orientation of incoming electron spin with respect to 
the Rashba electric field (Sec. llVjl . We summarize our 
main findings in Sec. 



II. HAMILTONIAN APPROACH TO 
QUANTUM SPIN-POLARIZED TRANSPORT 

A. Quantum coherence within the tensor product 
formalism 

In the general formalism of nonrelativistic quantum 
mechanics, spin degrees of freedom are internal, and are 
therefore described within a separate finite-dimensional 
vector space Tig which has to be multiplied tensorially 
with an orbital space Ho to get the full Hilbert space of 
quantum states 7i — Tio ® Ti-s- Any operator acting in 
Ti. can be expressed as a linear combination of the tensor 
products Oo Os, where Oo and Os act in Ho and Hs, 
respectively. A generic Hamiltonian model for spintron- 
ics H = Ho + Hs + Hso contains the following pieces: 
Ho = Oo ® which acts nontrivially in Ho and like a 
unit operator Ig in Hs] Hs = Io(E)Os acting nontrivially in 
Hs', and SO terms Hso which are linear combinations of a 
tensor product of spin and orbital operators coupling the 
two spin-polarized subsystems. The advantage of spin- 
tronics modeling through a quantum transport formal- 
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ism that takes as an input a microscopic single-particle 
Hamiltonian is that spin, as a genuine quantum attribute 
of an electron, is treated in the most natural way by an 
operator acting in a two-dimensional vector space Tig. 
The spin-dependent interactions are then accounted for 
by the terms in a Hamiltonian that are functions of the 
spin operator S = her/ 2 [a = {(Tx,CTy,(Tz) is a vector of 
the Pauli spin matrices] acting in TLs- An arbitrary state 
vector I^E*) G ?i is a linear combination of the separa- 
ble (i.e., uncorrelated) states forming a basis \4>a) ® jo"), 
where are basis states in TYq and corresponding ba- 
sis state |cr) in TLs are the spin eigenstates \a) of <t • u 
(ct =T,i is a quantum number labeling the two eigen- 
states, and u is a unit vector along the spin-quantization 
axis which is taken to be parallel to the magnetization of 
the ferromagnetic contacts in FM-Sm-FM geometry). 

Most of quantum-interference phenomena observed in 
standard mesoscopic experiments (such as, weak localiza- 
tion, conductance fluctuations, or the Aharonov-Bohm 
effects^) probe the orbital quantum coherence of electron 
states, which requires preservation of the relative phase 
of linear superpositions of spatial states of the electron.^ 
Such mesoscopic systems are smaller than the dephasing 
length (typically < 1 fim at low enough temper- 
atures T < IK) set by inelastic processes. Thus, co- 
herence times relevant for these investigations are com- 
pletely disconnected from the spin coherence times T2 
over which the relative phase in superpositions of f and 
I spin states is well-defined^ (transport of phase-coherent 
spin states has been observed^i over length scale L2 ~ 
100 ijm in semiconductors at low temperatures where T2 
can reach 100 ns). For an ensemble of spins created in 
experiments^-, the operational loss of phase coherence 
is described by the dephasing time T2* due to inhomo- 
geneities in the Zeeman terms (although such dephasing 
is independent of individual spin decoherence, it provides 
experimentally measurable lower bound T| < T2). The 
other basic time scale which, together with T2, captures 
essential features of spin dynamics is the spin relaxation 
time Ti (typically Ti > T2). The time Ti is classical in 
nature since it determines lifetime of an excited spin state 
(aligned long the external magnetic field), i.e., the relax- 
ation of unbalanced populations of spin states toward 
thermal equilibrium so that no concept of superpositions 
of quantum states is involved. ^■^ Corresponding spin dif- 
fusion length Li sets the upper limit to the size for any 
spintronic device since beyond this scale spin-encoded in- 
formation completely fades away^ The principal interac- 
tions providing modes of decoherence and spin relaxation 
are exchange coupling with other electron (or nuclear) 
spins and SO coupling to impurity atoms and defects. In 
spintronic systems of current interest, one can encounter 
quantum interference effects stemming from spin coher- 
ence2*24, both spin and orbital coherence that are inde- 
pendent of each other,'* or intertwined spin and orbital 
coherence when SO interaction arc strongiS^ 

These different quantum-coherent situations become 
transparent after being recast in the formalism of tensor 



products of Hilbert spaces and corresponding operators. 
In conventional macroscopic solids at finite temperatures 
one has to use a quantum statistical mechanics descrip- 
tion of electrons in terms of density operatorsSi (i.e., den- 
sity matrices when a particular representation is chosen) 

Pmacro = Po /Os. (1) 

which are determined by thermal equilibrium. The most 
general quantum state of a spin— ^ particle is accounted 
for by a density operatoriSLS^ 

Ps = l [is + p-o-y (2) 

where < |p| < 1. For example, |p| — 1 describes a pure 
quantum state which is fully polarized in the direction of 
vector p, while |p| = stands for a non-pure state (the 
so-called mixture) that is completely unpolarized. The 
trivial example is injection of both spin states in equal 
proportion, p, = 1/2(|T){T| + = t/^. Coupling 

of electron spin to the environment of surrounding elec- 
tronic or nuclear spins causes decoherence of spin quan- 
tum state, i.e., appearance of mixtures |p| < 1 of spin 
states instead of a single pure state. In semiclassical spin- 
tronics, the spin part of the full (non-pure) state becomes 
a pure state ps = Then, after taking only diag- 

onal elements of Po, a semiclassical description in terms 
of the Boltzmann distribution functions for f and J, elec- 
trons is obtained where phase-relations between differ- 
ent states are ignored^ (as applicable to systems with 
L > L^). Saying that the sample is phase- coherent in 
standard mesoscopic physics has a simple meaning in the 
formal language: an electron is described throughout the 
sample by a single wave function $(r) — (r|<I>) e Ho, 
i.e., the orbital factor gets "purified" po — |'I')('I'| at low 
enough temperatures where dephasing processes are sup- 
pressed below /im scales (while in usual experiments with 
unpolarized currents the spin factor is ps = 1/2 Is). In 
sample with both orbital and spin quantum coherence, 
the full state is pure p^^g^ — Pmcso, such as 

Pmeso= |$)(*|®|S)(S|, (3) 

The separability of the full state in Eq. Q holds only 
when scattering events are spin-independent (such as 
those generated by lattice imperfections, phonons, and 
non-magnetic impurities), thereby leaving the spin state 
of a traversing electron unchanged. In this case, the or- 
bital and spin factor states evolve independently in a co- 
herent fashion within its own factor subspaces Tig or Tio, 
under the unitary operator generated by the correspond- 
ing Hamiltonian. However, in the presence of spin-orbit 
interaction, these two quantum-coherences become inter- 
twined (note that if reduced density operator, obtained 
by tracing over the states in the other factor spaces, is 
pure, than it must be a factor of the total state opera- 
tor«2I). For example, if separable state with |t) spin factor 
is injected into a device with spin-orbit coupling, it will 
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evolve by unitary quantum evolution into a linear super- 
position I^E'i) of the separable state vectors \4>a) 'S) of 
some basis in H. Such superpositions are still pure, albeit 
non-separable quantum states, i.e., the different degrees 
of freedom (orbital and spin) of one and the the same 
particle have become entangled^ Thus, in this case in- 
dividual spin or orbital states lose their coherence since 
they are not pure any more (although no real ensemble of 
different quantum states corresponding to such mixtures 
exists in the samplej2S contrary to the case in the exam- 
ple of mixed spin states given above). This is analogous 
to the general decoherence mechanism that operates by 
entangling quantum system to its environment^ It also 
explains the formal meaning of 'spin-orbit quantum inter- 
ference effects' since entanglement is entailed by linear 
superpositions of states in the full Hilbert space H. More- 
over, such entangled states display quantum-mechanical 
non-locality, usually studied in the context of entangle- 
ment of two particles (e.g., Einstein-Podolsky- Rosen or 
Bell statesSi), which manifests here as "wholistic" spin 
effects all throughout the system^S When both 1 and J. 
particles are injected, the state of an electron is a statis- 
tical superposition |5'|)(\E'|| -|- rather than the 
coherent one -I- Thus, the quest for spintronic 
devices that would utilize quantum-coherent dynamics of 
all degrees of freedom requires to master the control of 
spins without destroying the phase coherence of orbital 
degrees of freedom^ (as is the case when many-body in- 
teractions with surrounding nuclear spins or other elec- 
trons^ are involved), and vice versam^ 

The effects of quantum-coherence on electronic trans- 
port are most easily discussed within the Landauer- 
Biittiker scattering approach. To study spin-polarized 
quantum transport, the basis of asymptotic scattering 
states at the Fermi energy Ep in the leads should be 
comprised of the following vectors 

{r\p)f = $p(y) ® exp{±ikpx) ® \a). (4) 

In the spinlcss case, the orbital factors are usually 
referred to as conducting ' channels '"iS; here a 'spin- 
polarized channel' pa has a real wave number kp > 0, 
a transverse wave function $p(r) defined by quantization 
of transverse momentum in the leads of finite cross sec- 
tion and given boundary conditions, and a spin factor 
state It) or ||). Thus, injecting a flux concentrated in 
the channel |in) = \p) (g) \a) will generate the following 
quantum state in the opposite lead 

lout) = ^tgp.^-„|g) (g) |cr')- (5) 

In the most general situation, with spin-dependent 
charge scattering present in the central region, this state 
is coherent mixture of all spin-polarized channels l@J. 
This equation introduces a transmission matrix t, whose 
elements are t^p^o-'o-- The t-matrix accounts for the uni- 
tary transformation that the central region would per- 
form on the incoming electronic wave function from the 
left lead. 



The formalism we introduce here is capable of deal- 
ing with all of the variety of quantum-coherent situa- 
tions relevant for transport (i.e., coherent superpositions 
of scattering quantum states in Ho, Ti-s, or Ho ^ Hs) 
in arbitrary geometry, and for arbitrary values of inter- 
action strengths, as long as one stays withing the do- 
main of single-particle picture. We choose the Anderson 
model for Ho (which is often used for numerical mod- 
eling of quantum transport^ localization-delocalization 
transitionjSi and brute force computation of conductance 
fluctuationsrSiS) 

Ho= ^em|m)(m|+ ^ t„„|m)(n| (6) 

\ ™ (m.n) / 

on a square N x N lattice (L = Na, with a be- 
ing the lattice spacing). Here tmn — te^'^"^'"", t is 
the nearest-neighbor hopping integral between s-orbitals 
(r|m) = ip{r — m) on adjacent atoms located at sites 
m = {jnx,rny) of the lattice. The local orbital basis is 
advantageous since it allows one to treat an arbitrary 
measurement geometry, which is an essential ingredi- 
ent in the studies of quantum transport through finite- 
size mesoscopic systems.^'* In the applied magnetic field, 
the link fields 0mn = — 0nm entering the Peierls phase 
factor introduce Aharonov-Bohm phases into the tight- 
binding model so that the flux through a given loop S is 
<I>S — X]{m n)eS '^mn in units of the flux quantum hje. 
The disorder is simulated by taking a random on-site po- 
tential with Em being uniformly distributed over the in- 
terval [—W/2,W/2] [the other possibilities to introduce 
disorder is random hopping, t £ [1 — 2W, 1], or random 
fluxes (pmn £ [0, Stt)]. The lattice site states of Ho multi- 
plied tensorially by |cr) define a basis |m) g) |cr) inH, and 
therefore a realCEDspin-space representation of the opera- 
tors in our algorithm. 

The purely spin part of H is exemplified by a Zeeman 
term 

Hs = -io ® • B, (7) 

where /i = g*^B/2 {g* is an effective gyromagnetic ra- 
tio and /is is the Bohr magneton) and B is the applied 
magnetic field. Finally, the relevant SO couplings in 2D 
heterostructures^ (neglecting oc terms) are given by 
the Rashba termi^ aR,m*(d-xv)-z/?i (z is a unit vector or- 
thogonal to the plane of the sample and m* is the effective 
mass), and the Dresselhaus term^ rjm*{vx(Jx — Vycfy)/h 
(generated by the lack of inversion symmetry of a peri- 
odic crystal potential in the bulk, like in GaAs): 

TT ctB.fi, ^ . . . ^ 

rjh ,^ - - - ^ 

where we replace m* by its tight-binding description in 
terms of t = /2m*a^ . The explicit use of the tensor- 
product notation in our formulas makes it possible to 
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obtain quickly and efficiently their matrix representation 
in a chosen basis. For example, to obtain the site-spin 
representation of the Rashba term in Eq. ((HJ we use the 
fact that velocity operator is a matrix 

{xn\vx\n) = -tmn{mx - n^) , (9) 

in the site representation, and then multiply, in the sense 
of direct product of matrices (when particular represen- 
tation is chosen, (X) stands for the Kronecker or direct 
product of matrices), this matrix with respective Pauli 
matrix. Thus, the final matrix elements in Eq. |(SJ) are 
just dimcnsionless constants multiplied by the material- 
specific "spin-orbit hopping parameters" = aR/2a 
and = ri/2a that set the energy scales of the Rashba 
and Dresselhaus terms, respectively. 

The Rashba term added to the tight-binding Hamil- 
tonian ((BJ will cause spin splitting [i.e., lifting of spin 
degeneracy for t and J, states having the same momen- 
tum, k = {kx,ky) ^ 0] of an energy subband of 2DEG 
into two branches 

) = Eo{kx,ky)±AE{kx,ky) 
= e + 2t[cos{kxa) + cos{kya)] 

±2tf^^Jsm'^{kxa) + sm^{kya), (10) 

for clean sample with constant e on the diago- 
nal [the splitting becomes linear^^ in the momentum 

AE = 2aY{,^kx + ky for parabolic subband dispersion 

Eo{kx,ky) — h{k^ + ky)/2'm* characterizing continuous 
systems] . Such splitting of the conduction band as a re- 
sult of spin-orbit coupling in the presence of an asymmet- 
ric confinement potential makes it a useful tool to model 
the electronic structure of confined narrow-gap semicon- 
ductors. Rashba SO term also induces, when viewed 
within the semiclassical picture, spin precession around 
effective Rashba magnetic field B/j(k) [Rashba Hamilto- 
nian, like any SO term, can be interpreted as the interac- 
tion of electron spin with k-dependent internal effective 
magnetic field, B(k) • S] with frequency w = AE/2h. 
In the presence of disorder, the precession around B(k) 
for a given k terminates after scattering of impurity (the 
other sources of scattering that change k are boundary 
surface and phonons). Then is starts again but along a 
different randomly selected axis, so that the change of 
spin direction by full precession, which occurs in ballis- 
tic systemSfSiii is suppressed by disorder. In the cases 
of our study where interplay of the Rashba interaction 
and charge scattering takes place (at zero temperature), 
initial full polarization of injected electrons is destroyed 
(technically, measurement of the spin-dependent proper- 
ties alone requires to use reduced density operator, which 
does not describe a pure state any more). This is analo- 
gous to D'yakonov-Perel' mechanism^ii which is usually 
discussed in the context of spin relaxation^'^, but in our 
study of transport through mesoscopic spintronics struc- 
tures, which are effectively at zero-temperature, it is re- 
sponsible for the spin decoherenceiS^ 



B. Landauer-type conductance formula for 
quantum spin-polarized transport 

The central result of our formalism is a direct algo- 
rithm to compute exactly the zero-temperature conduc- 
tance matri>i4 of a single sample 

G= (^^iT j =X^=xlvTi^ T^V' ^^^^ 

We imagine that the 2D sample is attached via per- 
fect Ohmic contacts (which are desirable but are cur- 
rently problematic to achieve at FM-Sm interface to- 
gether with efficient injection^^) to two ideal semi-infinite 
leads (Fig. ^) where one serves to inject spin-polarized 
current from ferromagnetic contacts (emitting t or J, elec- 
trons) into a semiconductor or nonmagnetic metal, and 
the other one drains electrons to the contacts detecting 1" 
and/or J, electrons. The leads also define, by transverse 
quantization, the orbital part of an asymptotic scatter- 
ing state in Eq. Moreover, they effectively remove 
"hot" electrons from the open sample, thereby bypass- 
ing the intricate issues of dissipation which must occur 
somewhere in the circuit to reach the steady-state regime 
under external pumping. While the trick of using ideal 
semi-infinite leads is a standard one in mesoscopic trans- 
port theory,™ the polarization of the asymptotic states 
stems from the different densities of states for f and J, 
electrons in the ferromagnetic contacts (which can be 
conventional metallic or semiconducting, currently avail- 
able only at low temperature^!) . The meaning of G is 
elucidated by casting it into a form analogous to the Lan- 
dauer formula^Siifl where the usual sum of the transmis- 
sion probabilities over all transverse propagating modes 
for spinless particles T = TrttMs replaced by a 2 x 2 ma- 
trix T of partial transmission coefficients describing the 
transition between left and right subsystems comprised 
of the two types ti i of spin-polarized asymptotic states 
Thus, TTT is the sum of transmission probabilities, over 
all conducting channels, for j" electron to travel from the 
left to the right lead whose magnetic moments are ori- 
ented in parallel. Analogously, T^J^ quantifies the proba- 
bility to detect (in the setup with antiparallel magnetic 
moments of the contacts) |-electron arising from spin pre- 
cessionSiii or spin-flips^ of the injected f-electrons. 

Diagonal components of G, G^^ and G-^-'-, are familiar 
from the studies of giant-magnetoresistancc'^* (in the case 
of spin-degenerate transport, trivially, G^^ = G^^ and 
G^-^ = G'^^ = 0). On the other hand, properties of the 
off-diagonal "mixing" conductances, G^^ and G'^'^ , are 
much less explored (also, for some problems of partially 
polarized transport between non-coUinear ferromagnets 
they have been introduced recently as complex quanti- 
ties^). The G-matrix allows one to compute the mea- 
sured conductance in different setups for spin-resolved 
experiments 

G"measurcd — • G ■ S^, (f ^) 
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where Si , Sc are vectors describing the type of incoherent 
mixture of f or J, electrons which are injected or collected, 
respectively. This quantum conductance expression can 
also be employed to get results of semiclassical spintron- 
ics by taking a disorder average and the semiclassical 
approximation in the large system limit, as has been the 
practice in related giant-magnetoresistance studies of the 
transport with two independent 'spin-channels'.^^ 

The conductance matrix Eq. l(TT|l is calculated by gen- 
eralizing a two-probe Landauer-type formula for the spin- 
degenerate transport, which is obtained in the linear- 
response and zero-temperature limit of an expression de- 



rived by the Keldysh technique 
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to a separate treat- 



ment of the two spin-polarized components 



N 

E 



(13a) 



t = 2A/-lmEi®/, -G^; 



N ■ 



-lmEi^ ®/s.(13b) 



The partial summation adds squared amplitudes of all 
elements of the transmission matrix t having the same 
spin indices 



t»i,TT 



t2(i-l) + l,20-l) + l, 



= t 



C2(i-l) + l,2j, 
2i,2(i-l) + lj 



C2i,2j, 



(14a) 
(14b) 
(14c) 
(14d) 



where i,j ~ 1,...,N. Here — ImS^jj ~ —{Y7^j^ — 
S2fl)/2i are nonnegative matrices with a well-defined 
matrix square root, where Sf^ ^ = [S£ are the self- 
energies (r-retarded, a-advanced) describing the "inter- 
action" of a sample with the left (L) or right (i?) lead. 
The 2N X 27V submatrix G^^ of the full Green func- 
tion matrix G'^nmaa' — (ill cr|G"'|m, cr') connects the lay- 
ers (i.e., rows of sites) 1 and N along the direction of 
transport (x-axis). The Green function, describing the 
propagation of t or J. electron between two arbitrary sites 
inside an open conductor in the absence of inelastic pro- 
cesses, is the site-spin representation of the Green oper- 
ator obtained by inverting the Hamiltonian for the given 
boundary conditions (hard wall in our case) 



(15) 



where Y7 = + ^/j- The self-energy matrices intro- 
duced by the leads are non-zero only on the end layers of 
the sample adjacent to the leads. Their analytical form 
is known exactly as S£^^(n,m) = i^g£_^(ns. His, tL), 
with 5£ ^(n5,m5,tL) the surface Green function of the 
bare semi-infinite lead between the sites rig and mg in 
the end atomic layer of the lead (adjacent to the corre- 
sponding sites n and m inside the conductor) iSiS^ Here 
and tc are the hopping parameters in the lead and on 
the lead-sample interface, respectively (see Fig. 



Although Eq. (O assumes dynamics of noninteracting 
quasiparticles, it is also valid for the zero-temperature 
linear-response regime of an interacting system^ (only 
single-electron elastic processes are allowed in this regime 
due to the requirements of energy conservation). In the 
scattering picture of transport of noninteracting quasi- 
particles on the tight-binding lattice, the orbital factor 
of the asymptotic states Eq. |0J is comprised of a quan- 
tized transverse wave function and a Bloch state (instead 
of a plane wave) , where the corresponding dispersion re- 
lation is defined by the semi-infinite tight-binding lattice 
modeling the leads. Since the Rashba electric field se- 
lects a preferred direction in space, one has to take into 
account the relative spin orientation of the incoming elec- 
tron with respect to the z-axis. We choose the represen- 
tation for the spin-part of the scattering state in Eq. Q 

as I T) ^ ) and I i) ^ ^ ^ 

rection of the spin-polarization of injected and collected 
electrons is defined by specifying the direction for which 
spin-operator is diagonal. For example, to study the in- 
jection of particles that are spin-polarized along cc-, y- 
and z- axis, we change the representation of the Pauli 
spin operators. For "f, J, along the z-axis 



Therefore, the di- 



a{z) = 



1 

1 



1 

-1 



(16) 



the usual textbook representation, where Gz is diagonal, 
is used; for spin-polarization in the y-direction 



a{Y) 



-i 

1 



1 

-1 



1 

1 



(17) 



and for spin-polarization along the x-axis 



1 
-1 



-i 

1 



-1 
-1 



(18) 



Note that only the diagonal matrix is unique, while the 
other two matrices are affected (as demonstrated here via 
one particular choice) by the freedom in choosing phase 
factors of eigenvectors forming the columns of unitary 
matrices U which transform the standard set of Pauli ma- 
trices into a new diagonal representation U^alJ (there- 
fore, it is more appropriate to talk about Pauli spin alge- 
bra whose operators satisfy abstract set of rules^i) . The 
injection of electrons which are spin-polarized in the di- 
rection of transport (cc-axis) and precess within the cen- 
tral region is a standard set-up of the spin-FET proposal. 
On the other hand, injecting y-axis (Fig^ polarized elec- 
trons, whose spin is conserved throughout the semicon- 
ductor, corresponds to spin valves^ 
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FIG. 2: Conductance Go (circles) and its variance Var Go 
(diamonds) in a time-reversal (B=0) and spin-rotation in- 
variant {tf^ = i° = 0) 2D conductor modeled on 32 x 32 or 
64 X 64 half-filled {Ef = 0) tight-binding lattice as a function 
of disorder strength W (disorder-averaging is performed over 
10000 samples). The dotted line is the value of the universal 
conductance fluctuations predicted perturbativelj"— for the 
metallic diffusive regime. 



III. CONDUCTANCE FLUCTUATIONS AND 
ANTILOCALIZATION OF SPIN-POLARIZED 
ELECTRONS IN DISORDERED RASHBA 
SPIN-SPLIT ELECTRON SYSTEM 



A major boost for the development of mesoscopic 
physics came from the pioneering experimental^ and 
theoretical^Sii^ studies of unexpectedly large conduc- 
tance fluctuations (CF), VarG ~ (e^/h)^ . Each phase- 
coherent sample is characterized by a fingerprint of re- 
producible (time-independent) conductance fluctuations 
as a function of magnetic field, Fermi energy, or change 
of impurity configuration at a fixed Ep. This violates the 
traditional notion that a given sample is well described by 
average values of physical quantities, particularly when 
Var G becomes of the same order as (G) [(...) denotes dis- 
order averaging] as L approaches the localization length 

Thus, even in metals, CF scale as Yar G/{Gf - L^^^d 
in c?-dimensions. Such counterintuitive behavior of CF 
and lack of self-averaging (especially in low-dimensional 
systems d < 2) necessitated the development of a "meso- 
scopic approach" as an alternative to standard statisti- 
cal mechanics description of macroscopic condensed mat- 
ter systems, where whole distribution functions of phys- 
ical quantities are to be studied in nanoscale quantum- 
coherent solids containing a macroscopically large num- 
ber of particles.'' In diffusive {£ ^ L metallic 
(G ^ e^/ft-) conductors, CF are considered to be univer- 
sal (UCF), VarG — Gd(2e^/ft,)^, where Cd is a constant 
independent of the sample size or the degree of disorder, 
at least within certain limits {Gd is reduced by a factor 
of 2, 4, or 8 by breaking time-reversal or spin-rotation 
invariance, or both^ for an exhaustive list of symme- 



try classes and possible crossovers in weakly disordered 
quantum dots and for weak SO coupling, see Ref. iASi) . 
However, numerical studies^^ of the evolution of Var G in 
3D mesoscopic metals show that CF monotonically de- 
cay with increasing disorder strength, becoming close to 
a constant predicted by the perturbative UCF theory^Sii^ 
only in a narrow range of W. Therefore, to have a refer- 
ence point for the subsequent study of fluctuations of G, 
we first calculate Var Go in the crossover between small 
and large W (i.e., from a quasiballistic, to a diffusive, 
and finally localized transport regime) for a conductor 
described by the non-interacting and spin-independent 
Hamiltonian Ho with B = = ^mn- The localization 
length in the band center Ep = of the paradigmatic 
Anderson model © is^i ^ ~ (1 -h 5.2 ■ 10'*/IF'*)a, while 
its mean free path at the same Fermi energy is given 
by ^ ~ 30a{t/W)^ . In this case, spin-degeneracy gives 
Go = GTT + G^^ = 2GTT (G^i = G^^ = q). The result 
plotted in Fig. 12 demonstrates that Var Go decreases sys- 
tematically with increasing W S [0, 10], while G2 = 0.186 
is the expected UCF value in 2D "metals" in the diffu- 
sive regime^ (in the quasiballistic regime CF are not 
expected to be universal^). 

A recent experimental investigation^ of transport 
through open ballistic (but chaotic due to the surface 
scattering) quantum dots in a GaAs heterostructure, that 
are exposed to a large in-plane magnetic field, point to 
the emergence of a Rashba SO interaction as the essen- 
tial ingredient in analyzing their CFm^t^ Inspired by 
their unexpected findings on the reduction of CF, we un- 
dertake a similar but spin-resolved analysis when spin- 
polarized electrons are injected into a 2D disordered open 
quantum dot. Thus, we "create" an experiment on the 
computer "measuring" the components of G for each re- 
alization of disorder, and thereby sample-to-samplc fluc- 
tuations VarGa/3 — {G\p) — (Ga/s)'^ as a function of the 
disorder strength W and Rashba spin-orbit hopping pa- 
rameter i^. Since the Rashba SO coupHng can be tuned 
in principle by an interface electric field we sweep 
(and neglect t^) from the week tf^ = 0.01 to the 
strong = 1-0 SO interaction limit^ (a similar range 
of Rashba hopping has been explored for spin-polarized 
transport through clean wires where it was shown that 
even in the ballistic case, subbands in Eq. ^| defined by 
the orbital factor states can be mixed depending on the 
strength of the Rashba coupling and wire widtbii). The 
important parameters to keep in mind are: tfo/W; the 
ratio £/L which delineates boundaries of different trans- 
port regimes in the crossover from ballistic £/L ^ 1 to 
diffusive £/L ^ 1 transport (e.g., the diffusive transport 
regime can be reached at much weaker disorder but in 
very large sample, while the interplay of coherent scat- 
tering off impurities and spin-dependent interactions is 
the pronounced when the energy scale of the terms in H 
introduced by Hg and Hso are comparable to the strength 
of the disorder); and £,/L since for ^/L < 1 electrons be- 
come localized. 

The mesoscopic fluctuating properties of the compo- 
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FIG. 3: Zero-temperature sample-to-sample fluctuations of 
the components of G for transport which is spin-polarized 
t, J, along: (a) the a;-axis, (b) the j/-axis, and (c) the z- 
axis. The 10000 samples are modeled on a 32 x 32 half-filled 
(Ef = 0) tight-binding lattice with diagonal disorder and a 
strong = 1-0 Rashba SO coupling. The circles on all three 
panels are conductance fiuctuations of the sum of the com- 
ponents of G, e.g., VarGz = VarfG^^ 4, qLI ^ 4, qLT^^ 
which describes injection and detection of both spin species 
polarized in the x-, y- or z-direction for VarGx, VarGy, and 
VarGz, respectively. 
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FIG. 4: Ratio Var Go /Var G of the conductance fluctuations 
in a pure disordered case (Fig. |5J to the fluctuations of the 
sums of the elements of G (G = Gx, Gy, and Gz) from 
Fig. 13 This reduction factor due to the Rashba interaction 
is plotted for: (a) the strong SO coupling limit tf^ — 1.0; and 
(b) the weak SO coupling limit tf"^ — 0.1 or 0.01. 



nents of G, for FM magnetization lying t, | in the x-, 
y- and z-direction, are shown in Fig. |31 for the largest 
tf^ = 1.0. The complete suppression of CF in the weakest 
disorder case (standard conductance fluctuations reach 
a maximum in the quasiballistic transport regime,'*^ as 
shown in Fig. is a feature of the strong SO coupling 
limit. While VarG^^^ ~ VarG-'"'' in all cases, fluctua- 
tions of the off-diagonal partial conductances Var G^^ 
and Var G^^ can display quite different patterns and even 
be independent of the disorder strength ("universality"), 
unlike the case of Go or the total conductances Gx , Gy , 
Gz introduced below. 

The counterpart of the measured conductance in ex- 
periments on unpolarized electrons'*^ would be to inject 

and collect both pure spin states Si = Sc = ^ | 

G = sJ ■G-Sc^^G^T + G^'^ + G^T-t-G^-^. (19) 
Therefore, Fig.^also plots the CF properties of the total 
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FIG. 5: The total conductance Gx,y,z = G" + G^^ + G^^ + 
G^^ describing injection and collection of both spin-species 
in the device of Fig. Q Disorder-averaging is performed over 
the same ensemble of samples studied in Fig. |3 Antilocal- 
ization effects Gx, Gy, Gz > Go are substantial (destroying 
Anderson insulator formed for > 6) at strong Rashba SO 
coupling tfa ~ 1-0, much smaller for = 0.1, and virtually 
gone at tf^ = 0.01. 



conductances Gx, Gy, and Gz- For example, G = Gz 
means that the spin part, |t) or |J,), of the asymptotic 
scattering state (@J is an eigenstate of ct^, i.e., fully po- 
larized along the z-axis. They should be contrasted with 
Var Go in Fig. |2] Comparison is facihtated by looking 
at the reduction factor Var Go /Var G in Fig. 0] which 
also depends on the direction of magnetization of the fer- 
romagnetic contacts (i.e., chosen spin-polarization axis). 
This is due to the breaking of rotational invariance by the 
Rashba term (whose electric field always lies along the z- 
axis). While Gx and Gz (characterizing the device hav- 
ing magnetization of the leads in the plane orthogonal to 
the direction of virtual Rashba magnetic field) are similar 
and exhibit a similar pattern of CF as a function of dis- 
order strength, the two setups can still be distinguished 
by looking at the CF patterns of partial conductances 
that are summed in Eq. (|19(l to get full conductances for 
unpolarized transport. 

Further inspection of the total conductances in Fig. [S] 
shows that below ly ~ 2, the system is in the quasibal- 
listic transport regime {i > L) where Go > Gx , Gy , Gz 
because of the additional scattering at the FM-Sm inter- 
face which can arise due to band structure mismatch^ 
induced by the Rashba term. Upon entering the diffu- 
sive regime antilocalization Go < Gz , Gy , Gz sets in as 
a standard quantum-interference effect (which survives 
disorder averaging) pertinent to all spin-orbit interac- 
tions affecting phase-coherent propagation through dis- 
order. ^4*^ Thus, strong SO coupling impedes localiza- 
tion effects on both the CF and conductances, even at 
very high disorder. At smaller SO couplings in Fig. |S1 
the antilocalization effects are substantially diminished 
(t^ = 0.1) or vanish Go « Gx,Gy,Gz altogether 
(*so = 0.01). Nevertheless, the presence of weak SO cou- 
pling (i.e., small tfo/W) is still palpable in the CF plotted 
in Fig. 21 and even more so in the relationship between 
different partial conductances studied in Sec. IIVI 



IV. SPIN-RESOLVED PARTIAL 
CONDUCTANCES IN DISORDERED RASHBA 
SPIN-SPLIT ELECTRON SYSTEM 

The long lifetime of the electron spin orientation has 
lead to plausible expectations that the conductance of 
a two-terminal device G^^ with parallel magnetizations 
in the ferromagnetic contacts should be higher than G^^ 
when that orientation is antiparallel^*Si However, a re- 
cent experiment^ has unraveled surprising result for the 
disorder-averaged conductance difference (G^^ — G^^): it 
decreases upon increasing sample length (InAs containing 
2DEG that is attached to permalloy leads), and eventu- 
ally turns negative in the quasiballistic transport regime 
I ^ L. Theoretical modeling^ of the effect of localized 
spin-flip interactions, which couple two spin subsystems 
at a discrete sets of points within a mesoscopic disor- 
dered sample, finds similar phenomenon under the sim- 
plified assumptions for quantum-mechanical transmissiv- 
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FIG. 6: Disorder-averaged relative conductance difference 
AG||/G|| as functions of the disorder strength (i.e., in differ- 
ent transport regimes: ballistic, quasiballistic, diffusive, and 
localized), Rashba SO coupling, and the direction of magneti- 
zation in the leads with respect to the z-axis Rashba electric 
field.. The injected electrons, into the same set of Rashba 
spin-split disordered conductors studied in Fig. |21 are spin- up 
polarized in the direction denoted on each panel. 



FIG. 7: Disorder-averaged relative conductance difference 
AG||/G|| as functions of the disorder strength, Rashba SO 
coupling, and the orientation of magnetization in the leads 
with respect to the z-axis Rashba electric field. The injected 
electrons, into the same set of Rashba spin-split disordered 
conductors studied in Fig. |21 are spin-down polarized in the 
direction denoted on each panel. 



ities (e.g., T^^) that are independent of the spin orienta- 
tion of the incoming electron. The results of Sec. lIIII point 
out that such a treatment is not sufficient in devices with 
a Rashba SO couphng (as hinted at in Ref.0). Therefore, 
we analyze in Fig. |3 the disorder-averaged conductance 
differences (normalized by respective sums) when an f- 
electron is injected 



(20) 



while Fig.Qplots the corresponding differences for injec- 
tion of a spin-down electron 



(21) 



using three different cases for the spin-polarizations of 
incoming electron with respect to the 2-axis Rashba elec- 
tric field. The set of disorder and SO hopping parame- 



ters is the same as in Sec. IIIII Comparing AG|j/G|| 

and AG||/G|| shows that the relationship between con- 
ductances for parallel and antiparallel magnetization in 
the leads is quite intricate: for the x-axis or z-axis mag- 
netization, the Rashba spin-split system is not invariant 
with respect to the spin-subsystem interchange; on the 
other hand, in the case of polarization along the y-axis 
invariance holds, AGj|/G|| « AG||/G||. We also find 
that some of the conductance differences for the x- and 
y-polarization direction can change sign, even for small 
SO coupling {tf^/W ~ 0.1). This occurs within differ- 
ent transport regimes, depending on the polarization of 
the incoming electron, so that negative values of conduc- 
tance differences are not confined only to the quasibal- 
listic transport regime as in Refs. I4ll2a (the most similar 
case here to their phenomenology is t polarized, along the 
z-axis, current injected into the conductor with moderate 
SO coupling tf^ = 0.1). 

These negative values for conductance differences in 
2D mesoscopic spintronics disordered systems is a mani- 
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festation of a novel quantum interference effect involving 
the two spin components.'* For instance, if an |-electron 
is injected into a ballistic sample, the Rashba interaction 
will induce spin precession rendering non-zero G^^ ^ 
because the |-electron could be detected at the drain even 
if no spin-down electrons are injected at the source i^iii 
When scattering off impurities is added to the SO inter- 
action, the injected electron pure quantum state, with 
definite spin and momentum, evolves into another pure 
quantum state that is a linear superposition (as discussed 
in Scc.^ of different and entangled momentum and spin 
eigenstates. In the intuitive (semiclassical) picture of 
Feynman pathsf2£ for some device configurations shown 
in Figs. El and |3 this superposition is equivalent to a 
destructive interference taking place along those trajec- 
tories where projection of the electron spin on the mag- 
netization axis of the source has not yet changed sign at 
the drain contact. 



V. CONCLUSION 

Using our efficient real(8>spin space formalism for the 
evaluation of partial spin-resolved conductances, ex- 
pressed in terms of a Landauer-type formula, we have 
investigated some of the standard interference effects in 
the quantum-coherent propagation of electrons through 
two-dimensional disordered conductor, which is attached 
to two ferromagnetic contacts and is exposed to a Rashba 
spin-orbit coupling. In the limit of a strong Rashba in- 
teraction, all computed conductances exhibit a critical 
value of ~ e^/h in a substantial interval of large dis- 
order, where Go for the time-reversal and spin-rotation 
invariant system is negligible because of strong local- 
ization effects. The freezing out of the fluctuations of 
the total conductances, VarGx, VarGy, and VarGz, is 
even more dramatic, in contrast to the pure disordered 
case where the CF decay exponentially fast. The im- 
pediment of the Anderson localization at strong disor- 
der resembles a weak-antilocalization effect due to SO 
scattering off impurities,*^ which is known to generate 
a metal-insulator transition in 2D systems (otherwise. 



they are Anderson insulators, albeit with an exponen- 
tially large ^ at weak disordes^i). Nevertheless, it has 
been known that the Rashba induced antilocalization can 
be considerably differen t ^°i^* from the standard weak- 
antilocalization expression of Ref. El While both con- 
ductance fluctuations and localization effects are thought 
to fall into three universality classes (determined only by 
the invariance properties of the Hamiltonian of a con- 
ductor with respect to time-reversal and spin-rotation 
symmetry operationsi^i), here we demonstrate how the 
tracking of individual spin subsystems (i.e., injection and 
detection of spin-polarized electrons) allows one to reveal 
non-universal features specific to the Rashba SO coupling 
effects on the phase-coherent propagation through disor- 
der. For example, one has to take into account the an- 
gle between the spin-polarization of an incoming electron 
and the Rashba electric field. Moreover, the relationship 
between spin-resolved partial conductances unearths ad- 
ditional quantum interference effects specific to the con- 
text of spintronics: for particular set of disorder, SO, 
and spin-polarization parameters, one can obtain spe- 
cial superpositions of the two spin components, which, 
being entangled with orbital states of mesoscopic disor- 
dered systems, generate higher conductance for the two- 
terminal device with antiparallel contact magnetization 
than for parallel configuration. 

Finally, wc emphasize that our findings are numeri- 
cally exact within the single-particle picture of transport, 
and therefore non-perturbative in both the strength of 
the disordesS: and SO interactionii. For example, at 
large enough W the concept of mean free path, which 
is relevant for semiclassical transport and perturbative 
quantum interference corrections to it, ceases to exisli^ 
{£ < a < \f)- Eventually, the localized phase is encoun- 
tered when L ^. 
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